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Abstract. Elliptic flow at RHIC is computed event-by-event with NeXSPhe- 
RIO. Reasonable agreement with experimental results on V2(tj) is obtained. 
Various effects are studied as well: reconstruction of impact parameter direc- 
tion, freeze out temperature, equation of state (with or without crossover), 
emission mecanism. 
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1. Motivation 

Hydrodynamics seems a correct tool to describe RHIC collisions however ^(f?) is 
not well reproduced as shown by Hirano et al. [U, ■ These anthors suggested 
that this might be due to lack of thermalization. Heinz and Kolb[ J2] presented a 
model with partial thermalization and obtained a reasonable agreement with data. 
The question addressed in this work is whether lack of thermalization is the only 
explaination for this disagreement between data and theory for V2(tj). 
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2. Brief description of NeXSPheRIO 

The tool we use is the hydrodynamical code called NeXSPheRIO. It is a junction 
of two codes. 

The SPheRIO code is used to compute the hydrodynamical evolution. It is 
based on Smoothed Particle Hydrodynamics, a method originally developped in 
astrophysics and adapted to relativistic heavy ion collisions [Oj- Its main advantage 
is that any geometry in the initial conditions can be incorporated. 

The NeXus code is used to compute the initial conditions T M „, j M and m m on 
a proper time hypersurface [ 4 . An example of initial condition for one event is 
shown in figure 1. 




Fig. 1. Example of initial energy density in the rj = plane. 

NeXSPheRIO is run many times, corresponding to many different events or 
initial conditions. In the end, an average over final results is performed. This 
mimicks experimental conditions. This is different from the canonical approach 
in hydrodynamics where initial conditions are adjusted to reproduce some selected 
data and are very smooth. 

This code has been used to study a range of problems concerning relativistic 
nuclear collisions: effect of fluctuating initial conditions on particle distributions [ 
|S] , energy dependence of the kaon effective temperature [0 , interferometry at RHIC 
[Ej, transverse mass ditributions at SPS for strange and non-strange particles [IE]- 
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3. Results 

3.1. Theoretical vs. experimental computation 

Theoretically, the impact parameter angle is known and varies in the range of 
the centrality window chosen. The elliptic flow can be computed easily through 

Jd 2 N/d^drjcos[2(cb-A)}d^ ^ 

< V2iv) >=< joJNim^ > (1) 

The average is performed over all events in the centrality bin. This is shown by the 
lowest solid curve in figure 2. 

Experimentally, the impact parameter angle ip2 is reconstructed and a correc- 
tion is applied to the elliptic flow computed with respect to this angle, to correct 
for the reaction plane resolution. For example in a Phobos-like way 

< v b / ec ( V ) >=< V ^ S{ll) > (2) 



^< cos[2(tA< - V 2 >0 )] 



> 



where 



and 



obs( s E l d 2 N/dcf )l dr ] co S [2( ( t> l -^ 2 )} 

v > {v) = Y^PN/dJdq (3) 

In the hit-based method, tp2° an d ^2° are determined for subevents r\ < and 
> respectively and if tj 2 is computed for a positive (negative) 77, the sum in tp2, 
equation 3, is over particles with 77 < (77 > 0). 

In the track-based method, t/i^ an d ^2° are determined for subevents 2.05 <| 
77 |< 3.2 and V2 is obtained for particles around < 77 < 1.8 and reflected (there is 
also an additional y/2 in the reaction plane correction in equation 2). 

In figure 2, we also show the results for v% bs (il) for both the hit-based (dashed 
line) and track-based (dotted line) methods. We see that both curves can lie above 
the theoretical < v\(rj) > (solid) curve. So dividing them by a cosine to get < 
v 2' VCC ( r l) > w iH make the disagreement worse: < v\(jf) > and < v^ Tec {rj) > are 
different. 

Since the standard way to include the correction for the reaction plane res- 
olution (equation 2) seems inapplicable, we need to understand why. When we 
look at the distribution d 2 N/d<frdri obtained with NeXSPhcRIO, it is not symmetric 
with respect to the reaction plane. This happens because the number of produced 
particles is finite. Therefore we must write 

^ = ^(77)[l + ^2^(r,)cos(n(0-0 b ))+^27;;f(77)sinK0-4))] (5) 
= v° bs (v)ll + E 2v °n bS ^) cosing - ^)) + £ 2v'° bs (r,) sin(77# - ^))](6) 




Fig. 2. Comparison of various ways of computing V2' solid line is using the known 
impact parameter angle fa, dashed and dotted lines is using the reconstructed 
impact parameter angle fa. 10PT stands for equation of state with first order 
transition, EbE, event-by-event calculation, FO, freeze out mechanism for particle 
emission. Data are from Phobos [0. For more details see text. 

It follows that 

v°2 bs (v) = v\{rf) cos[2(0 2 - fa)] + v' 2 b ( V ) sin[2(0 2 - fa)] (7) 

We see that due to the term in sine, we can indeed have < v% bs (r]) > larger than 
< 1*2(77) >, as in figure 2. (The sine term does not vanish upon averaging on events 
because if a choice such as equation 4 is done for -02, v'^iv) an d sin(2(^2 — <fib) 
have same sign. Rigorously, this sign condition is true if ip2 is computed for the 
same 77 as u 2 (??)• Due to the actual way of extracting ip2 experimentally, we expect 
this condition is satisfied for particles with small or moderate pseudorapidity.) In 
the standard approach, it is supposed that d N/dfair) is symmetric with respect 
to the reaction plane and there are no sine terms in the Fourier decomposition of 
d 2 N/d(f)dr) (equation 5); as a consequence, v 2 bs (?7) < v\(rj). 

Since the experimental results for elliptic flow are obtained assuming that 
d 2 N/d(f>dr] is symmetric around the reaction plane, we cannot expect perfect agree- 
ment of our < v\{rj) > with them. In the following we use the theoretical method, 
i.e. < v\{rj) >, to make further comparisons. 
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3.2. Study of various effects which can influence the shape of V2{rj) 

In all comparisons, the same set of initial conditions is used, scaled to reproduce 
dN/dr/ for T f . out = 135 MeV. 

First we study the effect of the freeze out temperature on the pseudo-rapidity 
and transverse momentum distributions as well as i>2 (jj) (this last quantity is shown 
in figure 3). We found that 112(77) and d 2 N/p t dp t favor Tt out — 135 MeV, so this 
temperature is used thereafter. 




Fig. 3. Comparison of ^2(77) f° r two freeze out temperatures. Abbreviations: see 
figure 2. 

We now compare results obtained for a quark matter equation of state with 
first order transition to hadronic matter and with a crossover (for details see 
We have checked that the 77 and pt distributions are not much affected. We expect 
larger w 2 for cross over because there is always acceleration and this is indeed what 
is seen in figure 4. 

We then compare results obtained for freeze out and continuous emission [111], 
Again, we have checked that the 77 and pt distributions are not much affected. We 
expect earlier emission, with less flow, at large \q\ regions, therefore, narrower 1*2(77) 
and this is indeed what is seen in figure 5. 

Finally, we note that compared to Hirano's pioneering work with smooth initial 
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conditions, the fact that we used event-by-event initial conditions seems crucial: we 
immediately avoid the two bump structure. To check this, it is interesting to study 
what we would get with smooth initial conditions. We obtained such conditions 
by averaging the initial conditions of 30 Nexus events. Again, we have checked 
that the 77 and p t distributions are not much affected but preleminary results shown 
in figure 6 indicate that now is very different, having a bumpy structure. The 
case of smooth initial conditions has a well defined asymmetry and the elliptic flow 
reflects this. The ellipict flow of the event-by-event case is an average over results 
obtained for randomly varying initial conditions, each with a different asymmetry. 
As a consequence, the average V2 has a smoother behavior but large fluctutations [ 
If0| and is smaller (around the initial energy density sharp peaks seen in figure 1, in 
each event, expansion is more symmetric. No such sharp peak exists for the average 
initial conditions). 




Fig. 4. Comparison of V2(rj) for first order transition (fOPT) and critical point 
(CP) equations of state. 
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4. Summary 

v 2 (ri) was computed with NeXSPhcRIO at RHIC energy. Event-by-event initial 
conditions seem important to get the right shape of V2(i]) at RHIC. Other features 
seem less important: freeze out temperature, equation of state (with or without 
crossover), emission mechanism. Finally, we have shown that the reconstruction 
of the impact parameter direction ip2, as given by eq. (4), gives vf" 1 ^) > v^f]), 
when taking into account the fact that the azimuthal particle distribution is not 
symmetric with respect to the reaction plane. 

Lack of thermalization is not necessary to reproduce v-i (r\) . The fact that there 
is thermalization outside mid-pseudorapidity is reasonable given that the (averaged) 
initial energy density is high there (figure not shown). A somewhat similar conclu- 
sion was obtained by Hirano at this conference, using color glass condensate initial 
conditions for a hydrodynamical code and emission through a cascade code [112). 



We acknowledge financial support by FAPESP (2004/10619-9, 2004/13309-0, 
2004/15560-2, CAPES/PROBRAL, CNPq, FAPERJ and PRONEX. 



5. References 
References 

1. T. Hirano, Phys. Rev. 065 (2001) 011901. T. Hirano and K. Tsuda, Phys. 
Rev. C 66 (2002) 054905. 

2. U. Heinz U and P.F. Kolb, J.Phys. G 30 (2004) S1229. 

3. C.E. Aguiar, T. Kodama, T. Osada and Y. Hama, J. Phys. G 27 (2001) 75. 

4. H.J. Drescher, F.M. Lu, S. Ostapchenko, T. Pierog and K. Werner, Phys. 
Rev. C65 (2002) 054902. Y. Hama, T. Kodama and O. Socolowski Jr., 
Braz.J.Phys. 35 (2005) 24. 

5. C.E. Aguiar, Y. Hama, T. Kodama and T. Osada, Nucl.Phys. A 698 (2002) 
639c. 

6. M. Gazdzicki, M.I. Gorenstein, F. Grassi, Y. Hama, T. Kodama and O. 
Socolowski Jr., Braz.J.Phys. 34 (2004) 322; Acta Phys. Pol. B 35 (2004) 179. 

7. O. Socolowski Jr., F. Grassi, Y. Hama and T. Kodama, Phys. Rev. Lett. 93 
(2004) 182301. 

8. F. Grassi, Y. Hama, T. Kodama and O. Socolowski Jr., J.Phys. G 30 (2005) 
S1041. 

9. B.B. Back, M.D. Baker, D.S. Barton et al., Phys. Re v. Lett. 89 (2002) 
222301. B.B. Back , M/D. Baker, M. Ballint ijn et al , |nucl-ex /0407012 

10. Y. Hama, R. Andrade, F. Grassi et al., |hep-p h/0510096| |hep-ph/051010'T| 

11. F. Grassi, Y. Hama and T. Kodama, Phys. Lett. B 355 (1995) 9; Z. Phys. C 
73 (1996) 153. 



12. T. Hirano, nucl-th/0510005 



8 



0.06 
0.05 
0.04 

u 

<n 0.03 
0.02 
0.01 




- 


i 1 i 1 i 1 i 1 

(15-25)% 

Average over 200 events - 




m m 1 i ,,: 






— 


s> \ 

\ \ _ 
■ \ \ 

\ \ _ 

■\\ 

\ \ 
\ ■ 

V* 


¥ 


■ exp. data 

— IOPTEoS (EbE-FO,T f = 135 MeV) 
-- lOPTEoS(EbE-CE) 




1,1,1 





Fig 



(15-25) % 




0.03 
0.02 
0.01 




■ Exp. data 






-- !OPTEoS(<IC> 


- FO, T f 


= 135 MeV) 


!OPTEoS(EbE- 


FO, T f = 


= 135 MeV) 




-2 





2 



-4 -2 2 4 

Fig. 6. Comparison of V2(rj) computed ev9nt-by-event (EbE) and with smooth 
initial conditions (< IC >). 



